Temperature-resonant cyclotron spectra in confined geometries 
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We consider a two-dimensional gas of colliding charged particles confined to finite size containers 
of various geometries and subjected to a uniform orthogonal magnetic field. The gas spectral den- 
sities are characterized by a broad peak at the cyclotron frequency. Unlike for infinitely extended 
gases, where the amplitude of the cyclotron peak grows linearly with temperature, here confine- 
ment causes such a peak to go through a maximum for an optimal temperature. In view of the 
fluctuation-dissipation theorem, the reported resonance effect has a direct counterpart in the electric 
susceptibility of the confined magnetized gas. 

PACS numbers: 05.40.-a, 52.20.-j, 76.20.+q 



I. INTRODUCTION 



Electronic fluctuations are of a great importance in 
plasma physics, due to their relevance in charge and en- 
ergy transport [l| and to the well-established connection 
between fluctuation spectra and electronic susceptibility 
0. As the temperature of the system is increased, the 
power of thermal fluctuations contained in any narrow 
frequency interval is also expected to increase. However, 
such a straightforward temperature dependence has been 
observed in systems where the electron dynamics is inco- 
herent, as a result of the electron interactions with other 
electrons, positively charged ions, and impurities. If, 
however, the electron dynamics also contains a coherent 
component such as, for instance, rotation with cyclotron 
frequency in the presence of a magnetic field, then the 
temperature dependence of the fluctuation power spec- 
tra may develop nontrivial resonant behaviors. For in- 
stance, it has been demonstrated Q that a magnetized 
plasma operated at the limit-cycle fixed-point bifurcation 
point and driven by a tunable external white noise un- 
dergoes, stochastic resonance ^^l. On the other hand, the 
observation that in constrained geometries the matching 
of thermal scales and characteristic system length can 
produce detectable resonant effects has been previously 
reported [5, 6]. 

In this paper we show that a simpler instance of tem- 
perature controlled resonance can naturally occur in a 
confined magnetized electron gas, due to the matching 
of two lengths, the electron intrinsic thermal length, or 
gyroradius, and the finite system size. The effect inves- 
tigated here should not be mistaken for a manifestation 
of the well know electron cyclotron resonance, which re- 
sults, instead, from the matching of two frequencies, the 
cyclotron frequency of an electron moving in a uniform 
magnetic field and the pump frequency of a perpendic- 
ular ac electric field Q- The dynamics of a magneto- 
plasma electron can be reduced to the two-dimensional 
(2D) Brownian motion of a charged particle subjected 
to a uniform magnetic field. In Sec. |lT] we analyze the 




FIG. f: (Color online) Power spectral density Sv^d'^) = 
Svyis^j), Eq. (|3]), for kT/m = 1, ljc = ^ and different 7 (re- 
ported in the legend). Inset shows the velocity correlation 
function of Eq. gj, G{t), for 7 = 0.1. 



power spectral density of a magnetized Brownian par- 
ticle moving in an unconstrained planar geometry. In 
particular, we notice that the amplitude of the cyclotron 
peak grows linearly with temperature. At variance with 
this remark, in Sec. IIIII our numerical simulations show 
that in constrained geometries the cyclotron peak goes 
through a maximum for an optimal temperature (Sec. 
nil Ap . which, in turn, is determined by the matching 
of system size and (temperature dependent) average cy- 
clotron radius (Sec. IIIIBp . In Sec. IIVI we also show 
that the observed resonant temperature dependence of 
the cyclotron peak is robust with respect to variations 
of the boundary conditions and the geometry of the sys- 
tem. Finally, in Sec. |V]we discuss possible applications 
of this effect to confined systems of magnetocharges in 
biological and artificial structures. 



II. UNBOUNDED ELECTRON GAS 

In an equilibrium neutral plasma electrons of charge 
q and mass m oscillate with characteristic angular fre- 
quency [1] Wp = [riQq^ /{m£Q)Y/'^ (in S.I. notation), where 
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no is the average electron density and ep the vacuum 
permittivity. In the following we restrict ourselves to 
weakly magnetized plasmas to ensure that the cyclotron 
frequency associated with Bq, Wc, is much smaller than 
Wp, i.e., ujc — qBo/iti <^ ujp. This allows us to reduce the 
dynamics of a magnetoplasma electron to the 2D Brow- 
nian motion of a charged particle subjected to a uniform 
magnetic field. This is a lon gsta nding problem in plasma 
and astroparticle physics |9l-[ll|. Here we limit ourselves 
to introduce the results relevant to the discussion of our 
simulation data. 

The corresponding Langevin equation reads 
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where the vector ^{t) — {^x{t),£,y{t)) represents two in- 
dependent Gaussian white noises with {Ci{t)) — and 

For numerical purposes, it is convenient to rescale both 
time, t — >■ ujct, and space, r — >■ r/A. We recall that ujc 
is the i3o-dependent cyclotron angular frequency and A 
is the T-dependent electron gyroradius A — \JkT j {muj'^) 
Q . In dimensionless units Eq. ([1]) reads 
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where bo is a unity vector parallel to Bo- Note that, 
in the absence of boundaries, the only free parameter 
in the dimensionless Langevin equation ^ is the scaled 
damping constant g = j/ujc- In the foregoing sections 
all results will be given in dimensional units for reader's 
convenience. 

For the linear and unconstrained dynamics of Eq. ^ 
the p.s.d. S{uj) = {\r{uj)\'^), with r(cj) standing for the 
Fourier transform of r{t), can be computed analytically 
by means of standard harmonic analysis [l^, • Taking 
advantage of the fact that (t) and {t) are uncorrelated 
white Gaussian noises, we obtain 
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where Si and Sy- denote the p.s.d. of the i = x,y compo- 
nents of the 2D vectors f and v, respectively. Note that 
at resonance lu — lUc, the peak of the transverse veloc- 
ity is Sy^iiJc) = v^J{2-f), with = = 2kT/m. 
In view of the discussion below, we recall here that 
the stationary autocorrelation function of the velocity, 
G{t-t') = {vx{t)vx{t')) = {vy{t)vy{t')), solely depends 
on the difference t—t' and is related to the inverse Fourier 
transform of the p.s.d. 5*^^ (w) through the Wiener- 
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For an unbounded planar electron gas, the station ary 
distribution density of the velocity is Maxwellian |9l-[ll| 
and does not depend on either the damping constant, 7, 
or the magnetic field. Bo, 
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The autocorrelation function of the electron coordinates, 
{r{t)r{t')), diverges as a function of t and t' because the 
free motion of electrons on the plane is unbounded. 

To this regard it should be noticed that the diffusion 
of a Brownian charge carrier on a plane perpendicular to 
a constant magnetic field is normal, that is, for asymp- 
totically large t, (r^) = 4Z3Bt with 
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This means that for > the particle diffusivity gets 
suppressed, as Db is smaller than the free diffusion coef- 
ficient Do = kT/{m-f). 

The typical velocity p.s.d., S'u^(w), and an example 
of velocity autocorrelation function, G{t), are depicted 
in Fig.[T] for different values of the damping constant 7. 
The peak at the cyclotron frequency ujc broadens as 7 
is increased. In the remaining sections of this paper we 
choose 7 to be much smaller than uJc, i.e. 7 <C This 
ensures that the cyclotron peak is well pronounced or, 
equivalently, that electrons with any given velocity v tend 
to perform many cyclotron orbits of radius Tc — v/uc, 
before being perturbed by the combined action of noise 
and friction. Accordingly, in the underdamped limit the 
diffusion coefficient Db tends to kTj/{muj'^) — A^7. 



III. FINITE SYSTEMS 

Our goal now is to compute the power spectral den- 
sities (p.s.d.), S'(w), of a confined 2D gas of electrons 
with finite temperature and for different geometries and 
boundary conditions. It should be noticed that S{uj) 
encodes important information about the electromag- 
netic transmission properties of the electron gas. In fact, 
S{lu) can be directly linked to the imaginary part of 
the gas electric susceptibility, via the fiuctuation- 

dissipation relation, S{uj) ~ kTlTa[K{uj)] 

We start with the simplest case of an electron gas 
trapped in a strip delimited by two walls parallel to the y 
axis and a fixed distance d apart. Electrons are assumed 
to be reflected elastically by the strip boundaries, which 
leads to a zero net transverse fiux in the x direction. 
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FIG. 2: (Color online) Transverse power spectral density, 
Sx{u}), of an electron moving in a strip of width d = 1 (see 
sketch in Fig. |3]), as computed from Eq. ([2} for 7 — 0.005, 
tJc = 1 and the scaled temperature kT/ra reported in the leg- 
end. Solid curves represent the corresponding p.s.d. for an 
unbounded planar system. Inset: amplitude of the cyclotron 
peak Sx{i->Jc) as a function of kT/m (solid curve) and its ana- 
lytical approximation, Eq. ([8} (dashed curve). Note that for 
= 1 the scaled temperature kT/m coincides with A^. 



We numerically integrated the dimensionless Langevin 
equation ^ for g in the range 10^"^ 10~^ and then 
restored dimensional units with Wc = 1 and fixed system 
size. With respect to the time unit, this corresponds to 
reporting uj in units of uJc, with no severe restriction on 
the actual value of Bq , but for the condition that Wc ^ Wp 
(see Sees. ITTl and IVl for more details). With respect to 
the space units, setting the width of the strip to a given 
value, d, makes the electron p.s.d. depend explicitly on 
the temperature, which had been eliminated from Eq. 
^ by expressing all lengths in units of A. Note that in 
our plots = 1, so that the scaled temperature kT/m 
boils down to the square of the gyroradius. A, introduced 
in Sec. HH 



A. Resonant cyclotron peak 

The typical p.s.d. of the transverse coordinate, Sx{i-o), 
are depicted in the main panel of Fig. [21 The finely dotted 
curves represent Sx{oj), as computed numerically through 
Eq. ([2|) with 7 = 0.005 and ideal refiecting boundaries lo- 
cated sX X — ±d/2 with d — 1. The solid curves are 
the corresponding Sx{!^) for an unbounded planar elec- 
tron gas, as predicted in Eq. (jS)). The height as well as 
the width of the cyclotron peak depend on the scaled 
temperature kT/m. 

Unlike for the case of an infinite system, the amplitude 
of the cyclotron peak, Sxi^^c), for the confined electron 
gas depends resonantly on the temperature, as shown in 
the inset of Fig.[2j As T is gradually increased, SxiijJc) 
goes through a maximum for an optimal temperature, Tc, 
whose dependence on the confinement geometry is inves- 
tigated in the following sections. This behavior may be 
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FIG. 3: (Color online) Trajectory samples of |«i(f), (a), and 
x{t), (b)-(c), for kT/m = 0.02, 7 = 10"^, d = 1 and a;^ = 1 
(black curves) or _Bo = (red curves). In (c) is a blow-up 
of the trajectory portion marked by a rectangle in (b) The 
horizontal dashed line in (a) represents {\v\) as computed from 
Eq. 0. 



reminiscent of stochastic resonance Q. However, we an- 
ticipate that here the optimal temperature is defined by 
the matching of two length scales, rather than two time 
scales, as it is the case in ordinary stochastic resonance 
[3, [l3l ■ Finite volume effects have been reported in the 
early stochastic resonance literature , but in a totally 
different context. 



B. Quantitative interpretation 

The resonant temperature dependence of the cyclotron 
peak can be qualitatively explained as follows. First 
we recall that, in the underdamped regime, the diffu- 
sion time of a charged Brownian particle across a strip of 
width d is strongly suppressed by the presence of a mag- 
netic field, especially for <^1. From Eq. (43) of Ref. 

i'^^) = hence, the transverse diffusion time in 

the strip can be safely assumed to be much shorter than 
the cyclotron period, as illustrated in panels (a) and (b) 

of Fig. m 

Then, we notice that, as the interaction of the electrons 
with the walls is elastic, the equilibrium distribution of 
their velocity is not affected by the boundary geometry 
and is still given by the Maxwell distribution of Eq. ([5|) . 

Moreover, the cyclotron radius of an electron moving 
with instantaneous velocity v \s Tc = v /lJc] in the regime 
of low damping, 7 <C Wc, we can neglect the effects of 
friction on the electron orbits. This means that the con- 
tribution to Sx from a circular orbit of constant radius 
is proportional to /2^ [see discussion following Eq. ^] , 
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that is, increases quadratically with v. However, this con- 
clusion apphes only to electronic trajectories with centers 
located a distance not smaller than away from the re- 
flecting boundaries. Indeed, when the electrons come too 
close to the boundaries, they repeatedly bounce off the 
walls, so that their trajectories get distorted [see Fig. 
[HJc)] . For an equilibrium distribution of the electron ve- 
locities, this surely happens when their orbit diameter is 
larger than half the strip width, 2rc > d/2. 

In view of the arguments above, we assume for simplic- 
ity that fast electrons with too large a cyclotron radius, 
say, Tc > d/4, do not contribute to the cyclotron peak, 
whereas only a fraction 1 — Arc/d of the slower electrons 
with Tc < d/4 do. On further noticing that from Eq. ^ 
the equilibrium distribution of v is 
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we obtain the following estimate for the amplitude of the 
cyclotron peak. 
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where T = AyJkT /ml [dujc) and Erf(. . . ) denotes the 
standard error function. Note that this estimate for 
Sx{<^c) systematically underestimates the corresponding 
simulation curve as we neglected the residual contribu- 
tion from electronic orbits larger, but not too larger, than 
d/A. By numerically evaluating Eq. ([5]), one concludes 
that both the resonance value of the cyclotron peak, 
<S'max(wc), and Tc grow quadratically with d, namely, 
kTjm ~ Gm{uJcdf and S'max(wc) ^ 2.8 • 10-3^2/7. 

We also stress that the emergence of a resonant cy- 
clotron peak is not conditioned by the elastic boundary 
assumption. Boundary randomness or fluctuations may, 
indeed, affect the residual contribution from large elec- 
tronic orbits with Vc > d/4, but not the bulk contribution 
estimated in Eq. 



IV. DEPENDENCE ON GEOMETRY AND 
BOUNDARY CONDITIONS 

Next we numerically compute Sx{iOc) as a function of 
the temperature for different geometries and boundary 
conditions, in order to demonstrate the robustness of the 
temperature resonance of the cyclotron peak. 

We consider here three conflning setups for the 2D elec- 
tron gas: 

(i) Box with internal semi-transparent wall. 
Consider a d x d square box, centered at the origin 
X — y = 0, and containing a semi-transparent internal 
wall, X — h with \h\ < d/2, parallel to the y axis. The 
internal wall acts like a porous filter letting charges pass 
through only if the x component of their velocity, Vx , is 
larger than a certain threshold velocity Vc, i.e. Vx > Vc- 
If Vx < Vc, the electrons are elastically reflected back 



into their half box. We Gx d = 12, h = and plot in 
Fig.lU^a) the amplitude of the cyclotron peak as a func- 
tion of the scaled temperature kT/m for different values 
of the threshold. For = the compartment wall is 
transparent, so that the effective width of the box is d, 
whereas for Vc > b the internal wall acts as an almost per- 
fectly reflective boundary, thus dividing the square box in 
two rectangular compartments of width d/2. Despite the 
different geometries, conflned cyclotron orbits are con- 
flrmed to generate a resonant temperature dependence 
of the electronic p.s.d., no matter what Vc- Moreover, in 
agreement with our analytical discussion of the infinite 
strip from Sec. IIIIBi the optimal temperature Tc, corre- 
sponding to the maxima of the plotted curves, diminishes 
with, — > d"^ /A:, by a factor 4 on increasing Vc from 
to oo and, thus, halving the container width. 

(ii) Infinite strip with semi-transparent walls. 

Let us consider the infinite strip of Sec. IIII Al with the 
important difference that now its parallel walls are semi- 
transparent, as described in (i). If Vx < fc the electrons 
are contained in the strip; if Vx > Vc the electrons exit 
one wall and reenter through the other one with peri- 
odic boundary conditions. In Fig.|31[b), the peak ampli- 
tude S{uJc) is plotted as a function of the temperature 
for different values of the threshold. Since the bound- 
aries are periodic for Vc — and reflective for Vc — oo, 
here our approximate estimate for S{ujc) from Sec. IIII Bl 
is expected to work well only as Uc — >■ oo. In this limit, 
Eq. ^ reproduces fairly closely both Tc and Smaxl^c). 
In the opposite limit of purely periodic boundary con- 
ditions, the cyclotron peak at resonance, SmaxC^c), gets 
enhanced, while Tc only weakly depends on Vc- 

(iii) Annulus with with a reflecting inner wall and (iv) 
with a semi-transparent inner wall. 

Next let the electrons be trapped in a circle with radius 
ri , which represent an ideal reflecting boundary. The in- 
ner space is divided by a second circle of radius r2, with 
T2 < fi, into two compartments. The inner circle is con- 
centric with the outer circle and its circumference works 
as a semi-transparent wall (case (iv)) with threshold ve- 
locity Vc (applied to the radial component of v). The 
case (iii) of an ideal reflecting inner circle, corresponds 
to setting the threshold velocity Vc — oo, see Fig.^c), 
the gas is confined to an annulus. As one can antic- 
ipate from the discussion in Sec. IIII Bl the maximum 
of the cyclotron peak, S'max(wc), decreases on increasing 
r2, see Fig.|4l[c). Correspondingly, the optimal tempera- 
ture, Tc, also decreases because the width of the annulus 
shrinks. The dependence of the resonant cyclotron ef- 
fect on Vc is illustrated in Fig.|Ud). The effect is most 
pronounced in the case of a perfectly transparent inner 
circle, Vc — 0. Indeed, lowering Vc enlarges the surface ac- 
cessible to the cyclotron orbits of the confined electrons. 
As suggested by Eq. ([5]), Tc and the maximum of S{ujc) 
grow quadratically with the effective transverse dimen- 
sions of the gas container; for the simulation parameters 
reported in Fig.^d), this corresponds to an increase of 
both quantities by a factor of about 4 as Vc raises from 
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FIG. 4: (Color online) Amplitude of the cyclotron peak, Sx{ijJc), as a function of the scaled temperature kT/m for the different 
geometries of Sec. IIVI as sketched. In all four panels Wc = 1 and 7 = 0.00166 (a), (b) and 0.003 (c), (d). (a) Square box 
containing a semi-transparent wall x = h [case (i)]. Simulation parameters are: d — 12, h = 0, and threshold velocity as 
reported next to the relevant data sets, (b) [case (ii)] Infinite strip with semi-transparent walls. Simulation parameters are: 
d = 18 and Vc as reported next to the relevant data sets, (c) Annulus with reflective inner boundary [case (iii)] with Vc — 00. 
The outer radius is kept constant, ri — 9, and the inner radius, r2, varied as in the legend, (d) [case (iv)] Same geometry as in 
(c) but with constant inner radius, r2 = 4.5, and semi-transparent inner wall with threshold velocity, Vc, in the legend. 



to 00. 

(v) Infinite strip with internal semi-transparent wall. 
Finally, we show that by appropriately choosing the ge- 
ometry of the system, the cyclotron peak can go through 
two maxima as a function of temperature. Such a dou- 
ble resonance was found by inserting an internal wall, 
X = h with \h\ < d/2, of tunable threshold Vc in the 
infinite strip of Sec. IIII Al Our simulation results for a 
symmetric geometry with d = 12, h — and different Vc 
are displayed in Fig. [5ja). The limiting regimes, Vc — 
and Vc ^ 00 are well reproduced by our approximate 
formula in Eq. ([S]). In the intermediate regimes, say at 
Vc = 2, the SicUc) features two maxima. The low temper- 
ature maximum is centered around the optimal tempera- 
ture, Tc°°\ corresponding to a partitioned strip, Vc = 00. 
Most remarkably, the optimal temperature of the second 
maximum on the right is systematically higher than the 
optimal temperature, Tc'^\ of the un-partitioned strip, 
Vc = ^- Moreover, such a double resonance could only be 
found for certain combinations of Vc and h/d, as shown 
in Fig.[5i;b). 

The occurrence of a double resonance can be explained 
by noticing that for electrons with v^ < Vc, the wall acts 
as an effective partition. In Fig.lSJa), where h = 0, this 
corresponds to splitting the strip into two equal strips of 



half width. Correspondingly, the slow electrons trapped 
in either half strip contribute a cyclotron peak that is 
the highest for T ~ Tc°°\ Fast electrons with Vx > Vc 
are free to move across the full width of the strip, so 
that the optimal temperature of their cyclotron orbits 
ought to approach Ti°' with ri°' ~ 4T^°°\ However, 
by taking a closer look at our derivation of Eq. ([8]), it 
is apparent that in the case of fast electrons the lower 
limit of the integral must be modified to account for the 
condition Vx > Vc- As a consequence of the bell shaped 
profile of the Maxwell distribution, such a modification 
of the integration range moves the optimal value Tc to 
appreciably higher values, in agreement with Fig.[SJ 



V. CONCLUDING REMARKS 

The temperature controlled resonance of the cyclotron 
spectra, emerging from a matching between the sys- 
tem size length d with the thermal electron gyroradius 
A, becomes detectable for a confined magnetized gas of 
charged particles under two important conditions, sum- 
marized by the inequalities 7 ^ Wc <C Wp. The condition 
Wc ^ Wp, introduced in Sec. |lTl required applying mag- 
netic fields of relatively low intensity. The underdamped 
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FIG. 5: (Color online) Amplitude of the cyclotron peak, 
Sx{u}c), as a function of the scaled temperature kT/m [case 
(v) of Sec. llVj for ojc = 1 and 7 = 0.00166. (a) Infinite strip 
of width d — 12 with a semi-transparent filter located at /i = 
and different thresholds Vc (in the legend). The curves of Eq. 
(|8]) for d = 12 (dashed) and 6 (dotted-dashed) are reported 
for a comparison; (b) Geometry as in (a), Vc = 2 and different 
h (in the legend). 



regime, 7 <i; Wc, was assumed to enhance the cyclotron 
peak of the transverse p.s.d., iS'x(w), over its background. 
Both conditions can be met in magnetoplasmas [ij. 

In normal metals the observation of the resonant cy- 
clotron effects reported here might seem out of question. 
At room temperature typical values of the electron damp- 



ing constant are 7 



10 
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or larger, so that an 



underdamped electron dynamics would set on only for 
exceedingly large magnetic field [l6j . A more promis- 
ing playground for an experimental demonstration of the 
effect under investigation is a 2D electron gas, where mo- 
bility can be quite high, thus, corresponding to a small 
damping constant, 7 ^ 10^ s^^. A relatively low mag- 
netic fields (of about 0.1 T) then would easily satisfy 
the condition 7 ^ ujc- However, since the plasma fre- 



quency, Wp, of an unconstrained 2D electron gas tend to 
be very low, artificial geometries should be implemented. 
To this regard it helps mentioning two sets of recent ex- 
periments, which detected, respectively, oscillations in 
the magnetoresistances of two-dimensional lateral surface 
superlattices with square patterns [l7j . and dynamical 
phase transitions between localized and superdiffusive (or 
ballistic) regimes for paramagnetic colloidal systems con- 
fined to magnetic bubble domains [lH. Both results can 
be explained, in semiclassical approximation, as a com- 
mensurability effect between the cyclotron radius of the 
magnetocharges and the spatial periodicity of the sub- 
strate, without the need to invoke quantum mechanics. 

We finally point out that the diffusion of confined 
magneto-charged particles is a topic of increasing inter- 
est not only in solid state physics. In medical research, 
for instance, magnetic nanostructures confined to 2D ge- 
ometries are thought to offer the most exciting avenues 
to nanobiomagnetic applications, including targeted drug 
delivery, bioseparation and cancer therapy, even if their 
diffusion properties are not yet fully controllable. The 
possibility of extending our analysis to nanobiomagnetic 
processes at the cellular level requires advances on at 
least two issues: (1) diffusion of complex magnetic ma- 
terials. In biomedical applications, pointlike magneto- 
charges are often replaced by synthetic magnetic struc- 
tures, such as magnetic microdiscs with a spin-vortex 
ground state [l^; (2) walls interactions. Contrary to 
our simple model, the interactions between a magneto- 
charge and cellular walls are typically inelastic, namely 
characterized by finite interaction times, energy trans- 
fer and even structural changes, like the activation on 
mechanosensitive ion channels. Both issues are the sub- 
ject of ongoing investigations by research teams world- 
wide. 
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